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ABSTRACT 

The extent to which the projected distribution of stars in a cluster is due to a large- 
scale radial gradient, and the extent to which it is due to fractal sub-structure, can 
be quantified - statistically - using the measure Q = fh/s. Here fh is the normalized 
mean edge length of its minimum spanning tree (i.e. the shortest network of edges 
connecting all stars in the cluster) and s is the correlation length (i.e. the normalized 
mean separation between all pairs of stars). 

We show how Q can be indirectly applied to grey-scale images by decomposing the 
image into a distribution of points from which fh and s can be calculated. This provides 
a powerful technique for comparing the distribution of dense gas in a molecular cloud 
with the distribution of the stars that condense out of it. We illustrate the application 
of this technique by comparing Q values from simulated clouds and star clusters. 



1 INTRODUCTION 



The dimensionless measure Q has been shown to be a ro- 
bust discriminator between clusters with a large-scale radial 
density gradient and clusters with small-scale subclustering 
(jCartwright fc Whitworth] |2004l |Schmeja fc Klessen||2006} 
|Elmeg rccn 20101. As it stands, the Q method can only be 
reliably applied to a collection of points, i.e. stars. However, 
given that star clusters are often embedded in gas clouds, it 
would be useful if Q could be adapted for grey-scale images, 
such as sub-millimetre maps. 

For a two-dimensional cluster of points, Q is equal to 
the normalized mean edge length of the minimum spanning 
tree (MST) fh divided by the normalized correlation length 
s. The value of s is defined as the mean separation between 
all points divided by the radius of the cluster. The MST is 
the shortest network of edges needed to connect together all 
the points in the cluster. The value of fh is its mean edge 
length normalised by the inverse square-root of the mean 
cluster surface density. Neither value by itself can distinguish 
between large-scale radial clustering and small-scale fractal 
sub-clustering; however, fh varies more with sub-clustering 
than s, and s varies more with radial clustering than fh. 
Because of this, the ratio Q — fh/s can distinguish between 
the two, with Q > 0.8 for radially clustered distributions 
and Q < 0.8 for fractally sub-clustered ones. 



Cartwright et al. ( 2006 1 have shown that the correlation 



length can be adapted for use on grey-scale images as the 
brightness of a pixel is analogous to the surface-density of a 
cluster. However, a robust grey-scale equivalent of the MST 
method has yet to be found. We present an alternative to 
directly analysing grey-scale images; this method involves 
decomposing the image into a collection of points, for which 
Q can then be calculated. 

In Section [2] we describe how a fh and s can be calcu- 
lated from grey-scale images. Section[3]details the models of 
star clusters and cloud images from which we calculate Q. In 



Figure 1. Starting from left to right, patches encompassing all 
pixels within n = 0, 1, 2, 3 and 4 pixel-widths from a chosen 
pixel. 



Section[4]we present and discuss Q values from both artifical 
clusters and clouds, conclusions are presented in Section [5] 



2 METHODOLOGY 

Suppose that we have a grey-scale image of N p i* = I X J 
square pixels, each with the same angular size SxS, and 
that the flux received from pixel is 



1 1.7 



Ida.. 



(2.1) 



pixel 



where I is the intensity (in whatever wavelength band is 
being used) and dfl is an element of solid angle. It follows 
that the total flux received is 



i=Xj=J 

i=i j=i 



(2.2) 



We decompose the image into A^ pnt points. The choice 
of N pnt is discussed later. The flux quantum is then chosen 
to be 

-ftot 



(2.3) 
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To convert the grey-scale image into an equivalent array of 
points, we pick a pixel lZ(i,j) at random, with no account 
taken of its flux Fij . 

If > AF, then we reduce 



TP. 



F^ — AF , 



(2.4) 



and place a point at Yij + Ar rn d, where r,j is the centre 
of pixel lZ(i,j) and Ar rll( j is a small random displacement 
(smaller in magnitude than the linear size of a pixel, S). 

If Fij < AF, we consider a patch of pixels. An n-patch 
comprises all the pixels whose centres lie at angular separa- 
tion < nS from the centre of pixel lZ(i,j); the configuration 
of n-patches for n = 0, 1, 2, 3 and 4 are illustrated on Fig- 
ure [T] We increase n until the flux from the n-patch exceeds 
or equals AF, i.e. 

F n _ patch = {Fij}>AF. (2.5) 

n-patch 

We then reduce the flux from each pixel within the n-patch, 
pro-rata, i.e. 

/ AF 

Fa -»■ Fij 1 



, ( 2 - 6 ) 

' n-patch t 

and place a point at position r pnt which is equal to the 
weighted centre of the removed flux, plus a small random 
displacement, i.e. 



Ipnt 



= E 



Fij Tij 



F 



n-patch 



+ Ar r 



(2.7) 



n-patch 

We repeat this process A^ pnt times, thus reducing F to t to 



zero: 



Ftot — > Fx 



=n p „ 

E 



Fo 
N, 



p nt 



(2.8) 



Note that for every iteration of this algorithm, lZ(i,j) is 
chosen completely at random and is thus permitted to have 
the same value more than once. Also, flux "detritus" left 



over from equations (2.4| and (2.61 is invariably swept up 



and accounted for by later iterations; the final iteration, for 
example, has an n-patch size which encompasses the entire 
image. 

From this collection of JV P nt points, we can now gener- 
ate the minimum spanning tree using Kruskal's algorithm 
( Kruskal|1956 1 with normalized mean edge length 

i=N-l 

where N = N pnt , An is the circular area of the point distri- 
bution (i.e. the smallest circle encompassing all points) and 
rrii is the length of edge i of the minimum spanning tree. 
The correlation length is then calculated 



N(N-1)R N 



=N-1 j=N 

E £i* 

i—l 3=4+1 



(2.10) 



where Rn = (Ajv/7t) 1//2 . Note that s is normalized by 
Rn and as the number of edeges within a two-dimensional 
area scales as iV 2 . For fh, the number of edges scales as 
N, thus the average edge-length needs to be d ivided by 
(NA) 1/2 (N - 1) (see |Cartwright fc Whitworth] \20Q4\ for 
more details). 



It can be argued that a circular area is not always the 
most intuitive shape to consider for these purposes; for ex- 
ample, some work with Q specifically considers non-circular 
distributions ( jCartwright fc Whitworth||2009[) an d different 
definitions of the area, ( Schmeja fc Klessen|20 06). However, 
as Rn = (Ajv/7r) 1//2 in equation (2.101, these area compo- 
nents cancel as Q = fh/s . 

By picking the intial positions of sampling cells ran- 
domly (rather than, for example, starting with the brightest 
pixel, as is done in some clump-finding algorithms) and fur- 
ther adding Ar rn( j to r pnt , we help to break up the lattice 
structure native to a grey-scale image and ensure that no 
two points lie directly atop one another. 

By averaging over patches of pixels, the algorithm 
smooths over small-scale intensity variations and hence lose 
information in converting images into points. These losses 
can be avoided by setting iV pn t to a sufficiently high value. 
For a grey-scale image with discrete integer values assigned 
to each pixel, setting iV pn t = F to t would enable an ex- 
act reconstruction of the image from its sampled distribu- 
tion of points. This would require approximately iV pn t ~ 
1 x J x 2 where iVbpp is the number of bits associ- 
ated to the value of each pixel. For a 100 x 100 pixel image, 
this would correspond to JV pn t ~ 10 6 for -/Vb pp = 8 and 
Npnt ~ 10 8 for A r b PP = 16 . Whilst this is valid for images 
with discrete integer value pixels, astronomical images are 
often composed of continuous floating-point value pixels. In 
this case, a distribution of discrete points can not exactly 
represent a grey-scale image for any practical value of -/V pn t. 

Through preliminary testing, we find that for an im- 
age generated with a flux-distribution function F(x,y), the 
resulting Q value is dependent on both image resolution 
and iV P nt, as shown in Figure [4] These dependencies can be 
largely mitigated by setting ratio of sampled points to pixels 
ATpnt/iVpix to a constant value, again indicated in Figure [4] 
Unless otherwise stated, we have set iVp n t 
results presented in Section [4] 



3 MODEL CLOUDS AND CLUSTERS 

3.1 Radial power- law distributions 

Centrally concentrated distributions of stars and interstellar 
gas can be constructed with with density 



P(r) = Pa 



(3.1) 



where p(r) is the density at radius r, po is a defined den- 
sity at fixed radius ro and a is the density exponent. A 
synthetic star cluster is created randomly using the Monte- 
Carlo method to position stars according to equation ( |3.1[ ). 
Such a cluster contains A* stars with positions 

r = ft r V(3-a) ; 

6 = cos" 1 (27l e - 1) , 

x = r sin(6 l ) co&{4>) , 
y = r sin(6 l ) sin(</>) , 
z = r cos(#) , 



(3.2) 
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where 7L r . IZg and 1Z& are random numbers between zero 



and one ( Cartwright fc Whitworth|[2004 1 . 

For a gas cloud with a density profile given by equation 
J3.ll, the surface-density at impact parameter b is 



E(6) = 2b 1 ' 



\6)dd. 







This can be calculated analytically for integer a: 

p(l-6 2 ) 1 / 2 ifa = 0, 
S(b)= i 2 ln 1+(1 - ft2)1/2 ifa = l, 

Lcos^b) if a = 2 . 



(3.3) 



(3.4) 



We can also solve equation (3.31 numerically for non-integer 



values of a. Using equation (3.4|, we can produce I x J 
pixel images by setting b — at the centre of the image and 
assigning each pixel a value of £(&). Examples are shown in 
Figure [3] 

3.2 Fractal distributions 

In contrast to centrally concentrated power-law clustering, 
multiscale sub-clustering in astronomy can be characterized 
using fractals. Fractals possess self-similar scaling defined 
by a fractal dimension D. Regular Euclidean geometry can 
also be shown to have a similar scaling property. For ex- 
ample, consider a cube of edge-length 1. If this is divided 
up into N — 8 sub-cubes then each will have edge-length 
I = 1/yN = 1/2. Assuming this relationship between N 
and I holds over all scales, it is said to have a fractal dimen- 
sion of 



D = 



log(AQ = Iog(8) 
log(l/i) log(2) 



= 3. 



(3.5) 



Now consider the same cube but this time populated with 
N — 4 sub-cubes of edge-length I — 1/2, i.e. four cubes 
and four cubic regions of empty space. If the sub-cubes 
are populated the same way over all scales, then it can 
be considered to have a fractal geometry with dimension 
D = log(4)/ log(2) = 2 . In general, any three-dimensional 
shape that scales with D < 3 is considered fractal (Voss 



1988). 



A synthetic fractal star cluster of dimension D can be 
constructed iteratively by considering a cube of edge-length 
2 with a parent-star at its centre. The cube is then divided 
into eight sub-cubes, of which a random 2 D are given a child 
star at their centre. The parent star is then deleted with the 
children-stars becoming parents such that the process can 
be repeated over N sell generations. A little noise is added 
to the final positions of the stars to break the cubic struc- 
ture and the distribution is pruned to a sphere of radius 1. 
This results in a fractally sub-clustered sphere of stars that 
is roughly self-similar down to a length-scale of 2 x 2~ JV|!on . 



Fractals are often used in this way to generate clusters ( Bate 



et al.|1998||Goodwin fc Whitw orth 2004| |Kouwenhoven et~ 



20101; conversely, as this paper details, there are also meth- 



ods for extracting fractal information from real observations. 

Following on from the synthesis of fractal star clusters, 
artificial fractal clouds of given D can also be contructed. 
We start by considering a cube of edge-length 2 and uniform 
density 1. For generation n = 1, the cube is split into 8 sub- 
cubes of I — 2 x 2~ n = 1, of which a random 2 D mature 




Figure 2. Generation of fractal cloud with D = 2.0 . Tiles show 
density fields projected onto the x-y plane. Begining with a uni- 
form medium in tile (a), we add fractal overdensities over six 
generations, as shown in tiles (b) to (g). Finallly, in tile (h) we 
crop the field into a sphere and project it on to an arbitrary plane. 



and are given a density p — 8 n = 8 such that p I 3 — 8 . The 
matured sub-cubes are recursively populated with sub-sub- 
cubes in the same way over N sen generations. In instances 
where 2 D is non-integer, the integer value is used to populate 
the current sub-cube and the remainder is passed on to the 
next sub-cube population. The process results in a cube di- 
vided into 8^™ sub-cubes of I = 2 x 2 _JVgon , approximating 
a fractal density field. The field is then pruned to a roughly 
spherical shape, i.e. p — for all sub-cubes outside a radius 
of 1. Note that throughout the recursive sub-division, the 
densities of unmatured sub-cubes are not set to zero. This 
serves to (i) introduce some residual background noise and 
(ii) avoid creating any "vacuum" within the density field. 

We then project the column density of these clouds 
onto a 2-dimension plane through a random line of sight 
to create an image. Image resolution is chosen to reflect 
the scale of self-similarity for chosen N SBn , i.e. the image 
size is 2 JVgcn x 2 Ngen pixels corresponding to 32 x 32, 64 x 
64, 128 x 128 and 256 x 256 pixels for 7V gcn = 5, 6, 7 and 8 
respectively. Figure [2] shows a step-by-step generation of a 
D — 2.0, six-generation fractal cloud. 



3.3 Perimeter-area method 

One of the most common methods for measuring the three- 
dimensional fractal dimension D from a grey-scale image is 
the perimeter-area method. This uses the relationship be- 
tween the perimeter P and area A of a two-dimensional 
fractal shape 

P oc A D2/2 , 



(3.6) 



where D2 is the fractal dimension in two-dimensions (Voss 



19881. This can be applied to iso-contour lines from a grey- 



scale image, where a plot of log(P) against log(A) produces 
a slope of D2/2. By measuring D2 of known distributions, 
D can be inferred (e.g. |Sanchez et al. (2005)). However, this 
method loses accuracy when D > 2.5 and can not algo- 
rithmically distinguish between large-scale central clustering 
and multi-scale sub-clustering ( Cartwright et al.|2006 l. 



3.4 Acquiring Q statistics 

We generate radial grey-scale images with 2.0 > a > 0.0 
and fractal grey-scale images with 3.0 > D > 2.0 for image 
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Figure 3. Point distributions and MSTs from radial and fractal grey-scale images. Row 1 shows the original grey-scale data, row 2 shows 
the distribution of points decomposed from the data and row 3 shows the MST connecting all the points. Column 1 has a radial density 
profile a = 1.5 , column 2 is uniform density and column 3 has a fractal dimension D = 2.5 . 



sizes of 32 x 32, 64 x 64, 128 x 128, and 256 x 256 pixels. 
Figure[4]illustrates how the dependency of Q on image size is 
minimized when we set N pn t/N p i^ — k where k is a constant. 
For simplicity and computational manageability we use k = 
1, as discussed in Section [2] 

The Q statistics from these images are shown in Table[l] 
and compared with those from artificial clusters. A hundred 
realisations are performed for five values of a and five values 
of D at each of the aforementioned image sizes. Examples of 
the grey-scale images, along with their sampled point distir- 



butions and minimum spanning trees, can be seen in Figure 

El 



The artificial star cluster Q statistics are generated from 
clusters of one hundred to one thousand points with specific 
radial and fractal distributions. As with the grey-scale im- 
ages, one hundred realisations are performed for each value 
of a and D. 



Statistical comparison of clouds and star clusters 5 



Table 1. Clustering statistics for radial-profile distributions. Note 
that N = 7Vpi x = JV pn t for all grey-scale images. Rows where 
N = "Cluster" pertain to artificial star clusters with one hundred 
to one thousand stars. Statistical uncertainties with apparant zero 
values are at least an order of magnitude smaller than those with 
finite values. 



Q vs Af pnt /JV pix 



1.00 


1024 


0.634 ± 0.005 


0.732 ± 0.001 


0.866 ± 0.007 




4096 


0.638 ±0.002 


0.732 ± 0.000 


0.871 ± 0.003 




16384 


0.639 ± 0.001 


0.733 ± 0.000 


0.873 ± 0.002 




65536 


0.638 ±0.001 


0.733 ± 0.000 


0.871 ± 0.001 




Cluster 


0.613 ±0.011 


0.735 ± 0.008 


0.833 ±0.015 


0.50 


1024 


0.646 ± 0.006 


0.774 ±0.001 


0.834 ± 0.008 




4096 


0.651 ± 0.003 


0.775 ± 0.000 


0.840 ± 0.004 




16384 


0.652 ±0.001 


0.776 ± 0.000 


0.841 ± 0.002 




65536 


0.650 ± 0.001 


0.776 ± 0.000 


0.838 ±0.001 




Cluster 


0.623 ±0.010 


0.775 ± 0.011 


0.804 ± 0.016 


0.00 


1024 


0.653 ±0.005 


0.807 ±0.001 


0.809 ± 0.006 




4096 


0.658 ± 0.002 


0.807 ±0.000 


0.815 ±0.003 




16384 


0.659 ± 0.001 


0.808 ± 0.000 


0.816 ±0.002 




65536 


0.657 ± 0.001 


0.808 ± 0.000 


0.815 ±0.001 




Cluster 


0.633 ±0.011 


0.808 ±0.008 


0.783 ±0.014 



Table 2. Clustering statistics for fractal distributions. 



D 



A' 



3.00 1024 


0.658 


± 


0.005 


0.809 


± 


0.001 


0.814 


± 


0.006 


4096 


0.660 


± 


0.003 


0.808 


± 


0.000 


0.817 


± 


0.003 


16384 


0.660 


± 


0.001 


0.808 


± 


0.000 


0.817 


± 


0.002 


65536 


0.658 


± 


0.001 


0.808 


± 


0.000 


0.815 


± 


0.001 


Cluster 


0.648 


± 


0.010 


0.809 


± 


0.010 


0.801 


± 


0.011 



2.75 1024 


0.615 


± 


0.020 


0.784 


± 


0.034 


0.785 


± 


0.023 


4096 


0.618 


± 


0.019 


0.787 


± 


0.034 


0.785 


± 


0.021 


16384 


0.612 


± 


0.021 


0.783 


± 


0.034 


0.782 


± 


0.022 


65536 


0.616 


± 


0.022 


0.792 


± 


0.039 


0.778 


± 


0.021 


Cluster 


0.609 


± 


0.022 


0.783 


± 


0.036 


0.779 


± 


0.025 



2.50 


1024 


0.574 ± 0.025 


0.769 ± 0.047 


0.748 ±0.032 




4096 


0.564 ± 0.028 


0.772 ± 0.054 


0.732 ±0.039 




16384 


0.553 ±0.024 


0.764 ± 0.050 


0.725 ±0.032 




65536 


0.551 ± 0.024 


0.767 ±0.058 


0.720 ±0.038 




Cluster 


0.570 ± 0.030 


0.777 ± 0.049 


0.735 ±0.043 


2.25 


1024 


0.516 ±0.030 


0.746 ± 0.072 


0.695 ± 0.045 




4096 


0.495 ±0.027 


0.744 ± 0.063 


0.669 ± 0.045 




16384 


0.471 ± 0.029 


0.740 ± 0.080 


0.641 ± 0.052 




65536 


0.464 ± 0.022 


0.751 ± 0.068 


0.622 ± 0.047 




Cluster 


0.494 ± 0.030 


0.737 ± 0.063 


0.673 ±0.052 


2.00 


1024 


0.478 ± 0.029 


0.739 ± 0.079 


0.651 ± 0.051 




4096 


0.443 ±0.027 


0.742 ± 0.080 


0.601 ± 0.049 




16384 


0.416 ±0.024 


0.737 ± 0.073 


0.569 ± 0.052 




65536 


0.394 ± 0.023 


0.744 ± 0.076 


0.535 ±0.054 




Cluster 


0.453 ±0.030 


0.739 ± 0.084 


0.619 ±0.067 



1.05 



1.00 



0.95 



a 


N 


m 


s 


Q 


0.90 


(b) 


2.00 


1024 


0.572 ± 0.006 


0.590 ± 0.001 


0.969 ± 0.009 








4096 
16384 


0.576 ±0.003 
0.577 ± 0.001 


0.586 ±0.000 
0.584 ± 0.000 


0.984 ± 0.005 
0.989 ± 0.002 


0.85 


(c) 




65536 


0.575 ±0.001 


0.583 ±0.000 


0.988 ± 0.001 








Cluster 


0.544 ± 0.012 


0.583 ±0.012 


0.932 ±0.020 


0.80 


(d) 


1.50 


1024 


0.611 ±0.005 


0.673 ± 0.001 


0.909 ± 0.008 






4096 


0.616 ±0.002 


0.673 ±0.000 


0.916 ±0.003 


0.75 






16384 


0.618 ±0.001 


0.673 ± 0.000 


0.918 ±0.002 








65536 
Cluster 


0.616 ±0.001 
0.589 ±0.012 


0.673 ± 0.000 
0.674 ± 0.008 


0.916 ±0.001 
0.873 ±0.016 


0.70 





0.05 



0.60 



0.55 



0.50 



0.45 



0.40 
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Figure 4. A plot of Q against the total number of points over 
the total number of pixels A/p n t/A/ p i x . The red squares represent 
Q values for 32 X 32 pixel images, the green circles for 64 X 64 
pixel images and the blue triangles for 128 X 128 pixel images. 
The group of lines labeled (a) represent Q values for a = 2.0 , (b) 
for a = 1.0 , (c) for a = 0.0 and D = 3.0 , (d) for D = 2.5 and (e) 
for D = 2.0 . Each data point is generated from fifty realisations. 



4 DISCUSSION 

Tables [l] and [2] list clustering paramemters for radial and 
fractal distributions. Column 1 lists the radial density expo- 
nent a and fractal dimension D. Column 2 lists the size of 
the images in pixels; rows were N = "Cluster" relate to arti- 
ficial star clusters. Columns 3 to 5 give the normalized mean 
edge length of the MST m, normalized correlation length s 
and Q = fh/s. Figures |6] and [7] show plots of m against s 
for 128 X 128 pixel images and artificial clusters respectively. 
As shown by Cartwright ( 2009]) this gives a second means 
of discriminating a and D by examining which area of the 
plot specific values of ffi and s fall into. 

Calculating Q does not provide direct information on 
the structure of a distribution. Instead, a and D are inferred 
by comparing a measured value of Q to that of a distribu- 
tion of known fractal or radial structure. An example of this 
is shown in Figure [5] where a fractal dimension or radial 
density exponent can be estimated by knowing the size of 
an image N and its measured Q value. 

The data in Figure [5] show how Q statistics from grey- 
scale images vary with image size. As a figure of merit, the 
closer the data is to dQ/dN — 0, the better the algorithm 
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Figure 5. A plot of Q against the total number of pixels N. 
Different data series represent different values of a and D. Note 
that the statistical uncertainties on the radial Q values compared 
with the fractal Q values are negligable due to the isotropy of 
radial profile distributions. 
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Figure 6. A plot of fh against s for artificial cloud data. Each 
data set represents one hundred realisations of 128 X 128 pixel 
images. 




Figure 7. A plot of m against s for artificial cluster data. Each 
data set represents one hundred realisations each with one hun- 
dred to one thousand points. 



scales with image size. We observe that for grey-scale data 
of 2.0 > a > 0.0 and 3.0 > D > 2.5 with image sizes of 
4096 > N > 65536 pixels, TV has no significant influence on 
estimating a or D. 

For images with TV > 4000 and 0.81 > Q > 0.73, we 
find that fractal dimensions in the range 3.0 > D > 2.5 
can be estimated from Q with an approximate uncertainty 
of <td ~ 0.1 irrespective of image size. For 0.73 > Q > 
0.54, estimating 2.5 > D > 2.0 carries a higher uncertainty 
of an ~ 0.2 and requires matching Q to TV as well as D. 
Radial density exponents in the range of 2.0 > ot > 0.0 
can be estimated when 0.99 > Q > 0.81 . We find that 
these estimates of a are largely independent of N, however, 
the statistical uncertainties are artifically reduced as, unlike 
the fractal distributions, the radial mass distributions are 
completely isotropic. 

When comparing Q statistics from both artificial clus- 
ters and images of artificial clouds, we often find that in dis- 
tributions of the same a or D, clouds have higher Q values 
than clusters. On closer inspection, it can be seen that this 
arises from higher values of m in the cloud data. This can 
be explained by considering the means by which point dis- 
tributions are constructed for clusters and clouds. For clus- 
ter generation, points are positioned using random numbers 



(see equation 3.2 \ and therfore are subject to Poisson "sub- 
clustering". When sampling points from images of clouds, 
whilst the position of the sample area is chosen at random, 
the sampling area extends over several pixels, thus averag- 
ing out some of the small-scale density variation. This goes 
some way to produce an anti-clustered distribution, which 
tends to lengthen the MST, thus increasing the value of m 
and Q. This can be seen in Figure |4j where an increase in 
TVpnt reduces the number of pixels from which a point is 
sampled, reducing anti-clustering and lowering Q. 

It is important to note, that these results relate purely 
to spherical artifical data. Elongation of star clusters has 
been shown to have systemattic effect on Q ( jCartwright fc| 
Whitworth||2009| |Bastian et al.||2009[ ), however this can be 
quantified and corrected for. It is reasonable to assume that 
this elongation effect also applies for greyscale data and will 
need to be considered in follow-up work. 
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5 CONCLUSIONS 

We demonstrate a method and analysis for taking Q mea- 
surements from grey-scale images of clouds. By decomposing 
an image into a distribution of points, we are able to apply 
the pre-existing methods of calculating Q and infer informa- 
tion on cloud structure. 

Whilst there are systematic differences between Q val- 
ues for clouds and clusters with the same a or D, this does 
not present a problem as the relation between Q and a or D 
can be calibrated independently for both types of data. We 
also find that grey-scale Q values are largely independent 
of image size for radial density profiles and fractal distri- 
butions with D > 2.5 . This makes Q a powerful tool for 
studying the structure of molecular clouds alongside that of 
star clusters. 

Q can also be applied to hydrodynamical simulations. 
By taking measurements at regular time-steps, the strutural 
evolution of both gas and sink-particle distribution can be 
quantified as a function of time. 
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